Towards better delineation of hydrothermal alterations via multi-sensor remote sensing and airborne geophysical data

Integrating various tools in targeting mineral deposits increases the chance of adequate detection and characterization of mineralization zones. Selecting a convenient dataset is a key for a precise geological and hydrothermal alteration mapping. Remote sensing and airborne geophysical data have proven their efficiency as tools for reliable mineral exploration. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Advanced land imager (ALI), Landsat 8 (L8), and Sentinel 2 data are widely-used data among various types of remote sensing images in resolving lithological and hydrothermal alteration mapping over the last two decades. ASTER is a well-established satellite in geological remote sensing with detailed Short-wave infrared (SWIR) range compared to visible and near-infrared region (VNIR) that controls iron-associated alteration detection. On contrary, ALI has excellent coverage of the VNIR area (6 bands), but does not possess the potentiality of ASTER for the SWIR and thermal regions. Landsat 8 is widely used and highly recommended for lithological and hydrothermal alteration mapping. The higher spatial (up to 10 m) resolution of Sentinel 2 MSI has preserved its role in producing accurate geological mapping. Notwithstanding the foregoing, implementing the four datasets in a single study is time-consuming. Thus, an important question when commencing an exploration project for hydrothermal alterations-related mineralization (orogenic mineral deposits in the current research) is: which dataset should be adopted to fulfill proper and adequate outputs? Here the four widely recommended datasets (ASTER, ALI, L8, and sentinel 2) have been tested by applying the widely-accepted techniques (false color combinations, band ratios, directed principal component analysis, and constrained energy minimization) for geological and hydrothermal alteration mapping of Gabal El Rukham-Gabal Mueilha district, Egypt. The study area is covered mainly by Neoproterozoic heterogeneous collection of ophiolitic components, island arc assemblage, intruded by enormous granitic rocks. Additionally, airborne magnetic and radiometric data were applied and compared with the remote sensing investigations for deciphering the structural and hydrothermal alteration patterns within the study area. The results demonstrated a different extent from one sensor to another, highlighting their varied efficacy in detecting hydrothermal alterations (mainly hydroxyl-bearing alterations and iron oxides). Moreover, the analysis of airborne magnetic and radiometric data showed hydrothermal alteration zones that are consistent with the detected alteration pattern. The coincidence between high magnetic anomalies, high values of the K/eTh ratio, and the resultant alterations confirm the real alteration anomalies. Over and above that, the remote sensing results and airborne geophysical indications were verified with fieldwork and petrographic investigations, and strongly recommend combining ASTER and Sentinel 2 results in further investigations. Based on the outputs of the current research, we expect better hydrothermal alteration delineation by adopting the current findings as they sharply narrow the zones to be further investigated via costly geophysical and geochemical methods in mineral exploration projects.

Irrespective of the utilized sensor and its capability, remote sensing data should be integrated with geophysical and/or geochemical data for better delineation of hydrothermal alterations-related mineralizations and result verification (through analyzing the harmony and matching among multisource datasets). Airborne geophysical data (e.g., magnetic, radiometric) has proven its efficiency and capability in detecting hydrothermal alteration tracts associated with mineral deposits. Aeromagnetic data allows delineating structural patterns that serve as channels for hydrothermal alteration and/or detecting altered zones connected to porphyry minerals [37][38][39] . Also, the aeroradiometric data via arithmetic combination (the ratio of two radioelement grids), for example, the ratio of potassium counts and equivalent thorium concentrations (K/eTh) or its inverted form 40,41 , the F-parameter (the ratio of the product of potassium (K) counts and equivalent uranium (eU) concentration to the equivalent of thorium (eTh) concentration) suggested by Efimov 42 , and Kd the deviation from the potassium (K) counts 43 are highly recommended for mapping potassic altered regions.
Consequently, one of the main targets of the current research is resolving the hydrothermal alteration patterns and structural framework of the study area based on multiscale (satellite remote sensing images, airborne geophysical data, fieldwork, and petrographic investigations) and multisource (multi-sensor remote sensing images with different spectral and spatial characteristics, magnetic, and radiometric grids) data to increase the probability of locating real alteration zones for further feasibility studies and exploitation. Additionally, the current research is an attempt to highlight a suitable dataset that well-matches ASTER findings in hydrothermal alteration mapping (as ASTER encounters some issues especially related to the availability of SWIR bands) through a comparative study performed with the widely utilized sensors for better identification of alteration zones. As using the four sensors to decipher the lithological and hydrothermal regimes of a particular study is time-consuming. We seek to define the true hydrothermal anomalies and assess their suitability and magnitude. Toward our aim, the well-known techniques, including false color combination (FCC), band ratio (BR), relative absorption band depth (RBD), directed principal component analysis (DPCA), and constrained energy minimization (CEM), are utilized. Additionally, our outputs are verified with aeromagnetic data (delineating the structural complexity), spectrometric gamma-ray data (highlighting the alteration zones), field observations, and petrographic investigations.

Study area and geological setting
The area of study includes Gabal El Rukham-Gabal Mueilha district, Eastern Desert, Egypt. They constitute the northwestern part of the Nubian Shield (Fig. 1a). This shield was developed through a prolonged series of orogenic events, referred to the Pan African Orogeny (~ 950 to 550 Ma 44,45 ), manifesting a superimposed tectonic framework of ductile and brittle deformations (D1, D2, and D3 46,47 ). This orogeny was accompanied by Neoproterozoic Pan-African thrusting (e.g., Beitan, Hafafit, Barramiya, El-Shalul, Meatiq), that is separated as suprastructural units from infrastructural gneisses, migmatites, and cleaved granites. Except for the southwest corner, covered by the Phanerozoic Nubian Sandstone, the majority of the study area is occupied by ophiolitic serpentinite, metagabbro, ophiolitic mélange, metavolcanics, and syn-to post-tectonic granite (Fig. 1b). This regime, which is characterized by diverse lithologies derived from various tectonic protoliths, is commonly associated with a number of mineral deposits (e.g., Au, Cu, and Sn) reflecting remarkable structural control and association with hydrothermal zones. The study area possesses various lithologies and favorable conditions for hydrothermal alterations and the associated mineralizations. In addition, it contains fifteen metallogenic spots 48 (Fig. 1b) that could be used as spatial verification points (besides the other datasets and field observations) for the results of hydrothermal alteration in the current study. Therefore, the area under consideration provides an excellent case study to evaluate the capability of ASTER, Landsat 8, Sentinel 2, and ALI for lithologic mapping utilizing the widely used FCCs and PCA in conjunction with the previous geological maps 24,49,50 .
Moreover, the ANS has structural and economic significance due to the presence of valuable base metals (e.g., gold, copper, zinc, silver, sulfides, chromite, nickel, etc. 44,45,[51][52][53][54] ). These precious metals, particularly gold deposits are linked with alteration zones in the central and southern Eastern Desert of Egypt (where the study area is located) developed by hydrothermal activities, induced by metamorphic or hydrothermal off-sets of unexposed intrusions (e.g., Barrmiya and Dungash 49,[55][56][57]. Several types of alteration zones are recorded in the  45 showing the location of the Barramiya-Mubarak belt (small black rectangle) where the study area is located, and (b) a geological map of the study area showing its lithological units, modified after Zoheir et al. 49  www.nature.com/scientificreports/ Barramiya-Dungash region, including listvenite, sericitization, sulphidization, carbonatization, and chloritization. In the Barrmiya area, gold mineralization is associated with listvenite, which developed along the sheared serpentinites and close to chromite, magnesite, and antimony ores 58,59 . While in the Dungash area, gold mineralization is associated with alteration zones of varying mineralogy superimposing a metamorphic assemblage of greenschist facies 56 .

Materials and methods
Remote sensing data characteristics and pre-processing. In this study, the ASTER, ALI, L8, and Sentinel 2 datasets are used to map lithological units and hydrothermal alterations. Table 1 shows the spectral and spatial characteristics of each sensor. ASTER was launched on the Terra platform in December 1999. In the current study, ASTER scene was obtained from the USGS (https:// earth explo rer. usgs. gov/). Earth Observing-1 (EO-1) started its mission on November 21, 2000, to provide multispectral data using ALI and hyperspectral data through Hyperion. For the current investigation, ALI data with a spatial resolution of 30 m for the multispectral bands and 10 m for the panchromatic band is applied. The Landsat 8 satellite was launched on February 11, 2013, and was equipped with two sensors, the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS), which allowed it to acquire high radiometric resolution (16-bit) data in the VNIR, SWIR, and TIR regions, as shown in Table 1. For the current purpose, a level 1 T (terrain corrected) scene was utilized. The Sentinel-2 satellite was developed by the ESA as part of the Copernicus Programme to provide spectral data in 13 bands with variable spatial resolutions (10, 20, or 60 m). For the current study, Sentinel 2A MSI, as an L1C product data, is accessed through the Copernicus Open Access Hub.
Unifying the data characteristics is recommended for such studies as the results may be affected by season (mainly due to the effect of rainfall), time of acquisition (that controls the shadows), and differences in pre-processing. In the current research, four cloud-free scenes were obtained beyond the range of flash flood events (After revising the timings of major flash flood events within the study area and its environs), with the following times of acquisition 07:59:31, 08:03:44, 08:13:29, and 8:16:11 for ALI (EO1A1740422003070110PZ), ASTER (AST_L1A_00303062007083043), Landsat OLI (LC81740432019298LGN00), and Sentinel 2 (S2A_ MSIL1C_20200505T081611), respectively. ASTER, ALI, and L8 are preprocessed (geometrically corrected according to UTM, WGS 84 zone 36 N, running Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes atmospheric correction) and clipped to the borders of the investigated area. Sentinel 2 data were georeferenced to the zone 36 North UTM projection using the WGS-84 datum, then atmospherically corrected using sen2cor processor in the Sentinel Application Platform (SNAP). All the datasets are resampled (nearest neighbor method) to a 20 m pixel size to allow comparison. A flowchart showing the adopted methodology utilized in the current study is introduced in Fig. 2.
Image processing methods. The current study utilized false-color composites (FCCs) and principal component analysis (PCA) for discriminating the different rock units and deciphering the alteration pattern within the study area. According to the data specifications, we focused on discriminating the geological units and comparing the results with each other through overlay analysis and visual comparison. Moreover, the results of the four datasets are compared to each other and correlated with a recently published geologic map of the  48 . Then band ratio (BR), relative absorption band depth (RBD), directed principal component analysis (DPCA), and constrained energy minimization (CEM) techniques have been applied to extract various types of alteration. These methods are adopted in the current research as they are widely accepted in the geological scientific community and their results are reliable according to several previous studies in different terrains 12,32,36,39,[60][61][62] .
Band ratio and relative absorption band depth. The Band ratio method is a powerful technique for discriminating different rock and mineral groups through boosting spectral variances and lowering topographical effects 63,64 . For iron-related alterations (sensed via the VNIR spectral range), ALI can detect all iron oxides by two different band ratios, while L8, ASTER and Sentinel 2 can detect them by a single ratio. The four datasets can map ferrous iron, ferric iron, ferric and ferrous oxides, laterite, gossan, ferrous silicates, as shown in Table 2.
To configure the whole alteration pattern within the study area, the well-known Sabins ratio 65 FCC is applied to the four datasets to compare the results with the known altered spots. In this ratio and mimicking Landsat TM results, hydroxyl-bearing alteration minerals are displayed in red, iron oxides (ferric iron in the case of ASTER) are shown in green, and ferrous oxides are represented by bluish colors. These anomalies are sliced and separated after specifying four threshold values related to each ratio statistic and according to the following equation: Th = μ + σ where μ is the mean value and σ represents the standard deviation. Due to its specified SWIR coverage (ASTER bands 5:7 are only covered by band 7 for L8, band 12 for Sentinel 2, and band 9 for ALI), ASTER data could specifically highlight areas invaded by other hydroxyl-bearing alterations (Table 2) compared to the general extraction of Sentinel 2 (11/12), L8 (6/7), and ALI (8/9). Consequently, relative absorption band depth (RBD) is applied to ASTER data for highlighting rocks enriched with alunite/ kaolinite/pyrophyllite ((4 + 6)/5), sericite/muscovite/illite/smectite ((5 + 7)/6), chlorite/epidote ((7 + 9)/8), and Si-OH ((6 + 8)/7) minerals 32 .

Constrained Energy Minimization method (CEM)
. CEM is a spectral mapping technique for the localization of minerals based mainly on employing a finite impulse response filter to allow only the desired targets to pass through while minimizing the output energy of any background noise according to a particular constraint 12 . CEM is a well-known approach and widely utilized with ASTER data 12,71-73 in geological and hydrothermal alteration mappings. To confirm the existence of alteration minerals and for more specifications, the CEM algorithm was applied to the potential ASTER SWIR bands 12,71-73 to verify and specify the dominant hydroxyl-bearing minerals (within the detected hydroxyl-altered zones). Utilizing laboratory reflectance spectra of muscovite, illite, kaolinite, pyrophllite, montmorillonite, epidote, and chlorite derived from the USGS spectral library, pixels enriched with these minerals are highlighted and visually compared with the previous results and the geological context of the study area. These minerals are selected as they are among the common hydrothermal alteration www.nature.com/scientificreports/ index minerals for propylitic alteration (e.g. chlorite, epidote), argillic alteration (e.g. kaolinite, illite, Montmorillonite), and phyllic alterations (e.g. Pyrophyllite). Additionally, all of them are OH-bearing minerals and could be used to check and verify the OH-bearing alteration zones in the current research.
Airborne geophysical data and processing techniques. The data used in this study is merged with the airborne magnetic and gamma-ray spectrometric survey performed by the Aero Service Division of Western Geophysical Company of America in 1983. The magnetic and radiometric data are available in the form of reduction to pole (RTP), potassium counts (K), and equivalent thorium (eTh) concentrations ( Fig. 4a-c, respectively). The data maps provided a clearer depiction of the boundaries between rock units and a more direct delineation of the magnetic and radiometric lineaments in the area under investigation, enabling for even more accurate identification of the types and positions of various structures (e.g., faults, contacts, and shear zones).
Since the main purposes of using the airborne geophysical data are to confirm the detected alteration pattern and map its structural framework, the following approaches are used to achieve this target: (a) Center for Exploration Targeting (CET) porphyry analysis, which is a new filter that contains a three stage sequential procedures (circular feature detection, boundary enhancement and tracing) to detect the porphyry magnetic signature of the hydrothermally altered rocks 37 ; (b) eTh/K ratio, which is derived from potassium and thorium concentration maps (Fig. 4b,c) to detect zones of potassic, phyllic and prophylitic hydrothermal alteration 40 ; and (c) 3D Euler deconvolution 74 , involves the orthogonal gradients of the potential field in the horizontal (x,y) and vertical (z) directions and can be expressed as where N is the structural index, B is the field base level, T is the magnetic field at x, y and z coordinates. The method applied on the RTP magnetic data (Fig. 4a) to map the structural framework of the altered areas and estimate their depths.

Results
Geological mapping. FCCs built from ASTER (9/4/1 in RGB) and Sentinel 2 (12/11/2 in RGB) data separate serpentinite rocks in a dark black color compared to greenish-yellow metavolcanics (shown in dark red color in the case of L8 (7/5/3 in RGB) and ALI (9/6/3 in RGB), as shown in Fig. 3. The selected bands for the FCCs are more or less around the same spectral range. For example, in the red channel of the given FCCs (displayed in Fig. 3), ASTER band 9, Sentinel 2 band 12, Landsat 8 band 7, and ALI band 9 are all representing SWIR of approximately 2.2 µm ( Table 1) that is very useful in geological discriminations. Similarly and after several trials, and checking previous studies related to the lithological mapping with these sensors 1,10 , the RGB composites for each sensor were selected. For lithological mapping, PCA is applied to all the bands of each sensor (shown in Table 1) to discriminate the lithological targets and compare them with previous geological maps and our field observations. ASTER PCs (2/1/3 in RGB) show serpentinites in a strong, golden yellow colour (confirmed in a  www.nature.com/scientificreports/ pinkish-white color in PCs 2/1/3 in RGB of L8) that could be easily delineated. Metagabbroic rocks are shown in a greenish colour for ASTER and L8 PCs. Unambiguous discrimination of the weathered syn-tectonic granites (greenish) from wadi deposits (greyish) is presented through Sentinel 2 PCs (2/1/3 in RGB) and confirmed through ALI PCs (3/2/1 in RGB) results (sky blue color for syn-tectonic granites). The sharp distinction between metavolcanics and serpentinites is also displayed through ALI and Sentinel 2 results. Separation of volcaniclastic metasediments (representing the mélange matrix) from the other rock units may be difficult as they contain spectral signatures of the neighboring rock units. Overall and after careful visual interpretation with previous geological maps and intensive field work, our results revealed that the four datasets are eligible (almost similar effectiveness) for highlighting all the exposed rock units (Fig. 3).
Hydrothermal alteration mapping. Results of the main hydrothermal alteration pattern (Table 3) revealed that the four-color composites of ASTER (4/7, 2/1, 4/3), Sentinel 2 (11/12, 4/2, 1/8), Landsat OLI (6/7, 4/2, 6/5),and ALI (8/9, 4/1, 8/6) are visually comparable, especially the hydroxyl-bearing alteration (Red channel), which exhibits almost the same pattern as shown in Fig. 5a-d. Compared to the geological maps, hydroxylbearing alterations are mainly confined to serpentinites. Metavolcanics are affected by iron oxides and hydroxylbearing minerals. Generally, iron oxides are mainly represented by blue and green colors and are restricted mainly to ophiolitic mélange represented by volcaniclastic metasediments, weathered syn-tectonic granites, and quaternary deposits. The known mineral deposits (Fig. 3a) are mostly located within the altered zones. This, in turn, confirms the eligibility of the four sensors for depicting the main hydrothermal alteration zones that are further verified and specified with other adopted techniques. It should be emphasized that relative band depths could highlight a wide range (category or group of minerals) of hydrothermal alteration and cannot specify a certain mineral when compared with spectral localization methods (e.g. CEM). For example, RBD could highlight the existence of phyllosilicates clay minerals (e.g. alunite/ kaolinite/pyrophyllite) 75 but cannot specify a certain mineral like CEM that could highlight a single mineral based on its spectral signature (reflectance or absorption spectra). Results of the RBD method applied to ASTER show the distribution of phyllosilicates clay minerals (alunite/kaolinite/pyrophyllite), Al/Fe-OH, Mg-Fe-OH, and Si-OH minerals as shown in Fig. 5e-h. A reasonable allocation of these minerals is notable concerning the exposed lithologies. For instance, alunite/kaolinite/pyrophyllite is mainly concentrated on intensively weathered granitic rocks in the northern part of the study area. Similarly, Al/Fe-OH-enriched minerals (mostly sericite/ muscovite/illite/smectite) are associated with syntectonic granite and ferruginous Nubian sandstone. Alternatively, chlorite/epidote minerals are approximately congruous with serpentinite rocks. Si-OH minerals are extensively represented over the ophiolitic mélange, metavolcanic, and serpentinites.
In Table 4, different kinds of alterations are shown based on how much the input band or band ratio contributed to the eigenvector matrix of DPCA. Scrutinizing the magnitude and sign of eigenvector loadings extracted from DPCA was extremely useful for delineating alteration minerals groups in the four datasets. For instance, in ASTER DPCA of bands 1, 3, 4, and 7, the first DPC does not reveal any specific information for any type Table 3. Thresholds, and anomalies of the detected hydrothermal alterations (and their intersections), using band ratios of the four utilized datasets.  (Table 4). On contrary, DPC 4 has strong positive loadings for B7 (0.96) and very small positive or negative loading for bands 1, 3, and 4. Thus, DPC4 (Fig. 6a) strongly reveals the desired target from this combination (hydroxyl-bearing minerals). For Landsat 8, DPC 1 is less informative due to the dominance of positive loadings for all the input bands. DPC 3 (Fig. 6b) is better than DPC 2 and DPC 4 in distinguishing hydroxyl-bearing minerals due to its strong loadings for band 7 (0.67) compared to the other components (Table 4). In the Sentinel 2 findings, DPC1 doesn't introduce any unique information due to the dominance of negative loadings for all the input bands. DPC 2 has moderate loadings of approximately all alteration mineral groups and thus could not strongly highlight any of them. DPC 3 could reasonably highlight iron oxides due to the negative loadings in band 2 and the positive ones in band 8a; however, discriminating hydroxyl-bearing minerals is extremely difficult due to the   (Fig. 6c) could distinguish hydroxyl-bearing minerals due to strong loadings in bands 11 and 12. Similarly, ALI DPC 3 (Fig. 6d) is the best to highlight OH-enriched minerals due to the high loadings in band 9 compared to the same band loadings in the other components. Consequently, our results are greatly harmonized with the previous studies 7,76-78 that accentuated the presence of the altered anomalies in the noisiest components (the third or the fourth), and the current research revealed that the best DPCs for disclosing hydroxyl-bearing minerals are ASTER DPC 4, Sentinel 2 DPC 4, L8 DPC 3, and ALI DPC 3. In the same way, by analyzing the eigenvector matrix to distinguish iron oxides' bearing rocks, ASTER DPC 2 (Fig. 6e), L8 DPC 3 (Fig. 6f), Sentinel 2 DPC 3 (Fig. 6g), and ALI DPC 3 (Fig. 6h) are the most informative components. For Sentinel 2, L8, and ALI, the desired wavelength (b2, approximately 0.48 µm) for distinguishing iron oxides has a strong loading in DPC 3; despite the lack of this wavelength in ASTER, DPC 2 was the best component, exhibiting the strongest loading (0.71) in band 1. DPCA was also applied to informative band ratio indices to specify the alteration types 61 . For ASTER data, the inputs were 2/1, 4/2, ((5 + 7)/6), and ((7 + 9)/8) to highlight rocks enriched with ferric/ferrous oxides, Al/Fe-OH, and Mg-Fe-OH. In DPC 3, ferrous oxides' bearing rocks (Fig. 7a) could be mapped due to the strong contribution of the B2/B1 ratio input (0.99) compared to the other inputs. DPC 4 discriminates lithological units enriched with ferric iron (Fig. 7b) due to the strong negative loadings in B4/B2. Our results revealed that DPC 1 (Fig. 7c) is extremely useful in discriminating Al/Fe-OH rich minerals (5 + 7)/6) from Mg-Fe-OH enriched minerals ((7 + 9)/8). In this way, ASTER can not only detect hydroxyl-bearing minerals but also classify them as Al/Fe-OH or Mg-Fe-OH minerals, reflecting the absolute superiority of ASTER data due to its SWIR bands. Moreover, using the CEM technique, some of these minerals could be nominated and specified. Results of Sentinel 2 DPCA revealed the sublimity of DPC 3 (Fig. 7d) in recognizing rocks enriched with OH-bearing minerals due to the strong contribution of the 11/12 band ratio (0.98). The differentiation between ferrous (11/4) and ferric (11/8a) rocks was clear in DPC 4 (Fig. 7e), where the first has reasonable negative loadings (− 0.52) compared to the strong positive loadings (0.82) of the second. Analyzing the DPCA results of L8 data revealed the superiority of DPC 3 (Fig. 7f) in the distinction of rocks enriched with OH-bearing minerals due to the strong contribution of the 6/7 band ratio (0.99). The distinction between ferrous (6/4) and ferric (6/5) was given through DPC 4 loadings (Fig. 7g). Similarly, ALI discriminates rocks enriched with OH-bearing minerals in DPC 3 and collectively highlights iron oxides in DPC 4 (Fig. 7h).

BR
As a continuation of revealing the hydrothermal alteration patterns in the study area, the CEM technique is implemented in the current study to approximately specify alteration minerals in SWIR (using ASTER data) ranges. CEM results and as an additional validation method (Fig. 8) showed that approximately all the highlighted pixels for muscovite, illite, kaolinite, pyrophyllite, montmorillonite, epidote, and chlorite minerals are confined to the highlighted hydroxyl bearing zones (previously detected using DPCA and BRs). Moreover, chlorite and epidote are concentrated along serpentinite rocks. Clay minerals (e.g., montmorillonite, illite, and kaolinite) and muscovite are mostly restricted to weathered granitic rocks or wadi deposits. These results not only confirm the highlighted zones and the efficiency of the utilized techniques but also highlight ASTER data for specific alteration mapping. Therefore, results from the other three sensors were compared to ASTER using anomaly overlay analysis to select the best sensor to be combined with ASTER for further geological mapping. Integrating with ASTER data gives confirmation for OH-bearing altered zones (mostly associated with hydrothermal mineral deposits), where the other sensors have different bandwidths and wavelengths. However, before comparing the areal extents of all the detected hydrothermal alteration anomalies, a comprehensive airborne geophysical analysis was applied to evaluate the reality of these detected zones using airborne geophysical data, including the reduction to pole (RTP) map, potassium concentration (K) map, and equivalent thorium concentration (eTh) map.
By collaborating with magnetic and radiometric data, which are commonly used to map magnetic and radiometric features associated with hydrothermal alteration rocks. We could confirm the recording of shear zones and other related structures that more or less coincide with the formerly mentioned hydrothermal alterations in all interpreted maps of the investigated area, besides estimating their depths. Concerning aeromagnetic data, a set of algorithms for CET porphyry analysis are applied to the RTP magnetic map to locate the porphyry-like circular or semi-circular structures (intrusions) that are associated with hydrothermal alteration. The intrusions and/or the inner alteration zone are associated with high magnetic anomalies, while the outer ones are less magnetic. In the preliminary phase, the circular features and their centers in the study area are highlighted by applying the Circular Feature Transform (CFT) (Fig. 9a), and then the Amplitude Contrast Transform (ACT) is used to delineate the boundaries of the detected circular intrusions (Fig. 9b). The yields of the CET and ACT  (Fig. 9c). On the other hand, the aeroradiometric data in the form of equivalent potassium and thorium concentrations assist to discriminate between different rock types in the study area; for example, ophiolitic serpentine is characterized by low values in K and eTh maps while the highest values are related to syn-orogenic granitoids. In addition, the highest values of the K/eTh ratio map (Fig. 9d) characterized by pink color, reveal the areas of hydrothermal alterations. The visual inspections of Fig. 9a-d reveal the areas of much being hydrothermally altered, which are associated with ophiolitic serpentine www.nature.com/scientificreports/ and oriented in an ENE-WSW direction over Gabal Um Salatit and Gabal Um Salim (the Baramiya-Mubarak shear zone), which is well-explained in the following section.
The structural context of the detected anomalies. Fowler 79 hypothesized that structural discontinuities in the Central Eastern Desert, such as major transpressional and transtensional shear zones and moderately dipping thrusts, aided in the intrusion of syn-tectonic granitoid intrusions, many of which are linked with hydrothermal activity and mineralization [80][81][82] . Among these structures, the northwest transpressional Najd-style and its conjugate Barramiya-Mubarak shear zone of the spectacular northeast trend and their second-and thirdorder extensional shears, which are genetically linked with ore deposits, particularly gold mineralization 29,81,[83][84][85][86] . Gold mineralization is claimed in a variety of distinct rock types, including ophiolitic serpentinite, arc metavolcanics, gabbroic, and granitoid rocks 87 . The Baramiya-Mubarak shear zone is detected locally in the research region at Gabal Um Salatit, which is composed of ophiolitic rocks and has a dextral displacement in the ENE-WSW direction (Fig. 10). Although this shear zone and its secondary fractures are not directly responsible for mineralization, they provide a pathway for the transformation of hydrothermal fluids, which are mainly responsible for the detected hydrothermal features in the current research.
To better assist in mapping structural features (key elements for hydrothermal alterations) of the investigated area and estimate their depths, the aeromagnetic RTP map (Fig. 4a) is used to perform the 3D Euler deconvolution method. The method was utilized with structural indexes 0, 0.5, and 1 for contacts, dykes, and faults, respectively. It is observed that the main tectonic trends in the area are oriented in NW-SE and ENE-WSW directions. The obtained depths are illustrated in Fig. 11a-c. A statistical analysis was performed to determine the prevailing depths. It is shown that most of the depths fall in the range of 0-600 m, indicating shallow sources, while the deeper depths up to 1500 m indicate moderate and deep depths of dykes and fault-like structures. In conclusion, the aeromagnetic data in the form of RTP combined with the Euler deconvolution clearly enhanced the appearance of the shear zone and other magnetic structures (contacts, dykes, and faults) that mostly constitute the pathways for hydrothermal fluids and gave an estimation for their depths.

Field verification and petrographic investigations
The main finding of the current research is highlighting the higher capability of collaborative usage of ASTER and Sentinel 2 data for better delineation of hydrothermal alteration zones. Our results coincided with airborne magnetic and radiometric data investigations, known mineralized zones, and previous studies. Additionally, fieldwork was performed to verify the detected hydrothermal alteration zones detected by the four datasets (Figs. 12, 13, and 14), and displayed in Fig. 15a to manifest the power of blended ASTER and Sentinel 2 (Fig. 15b) in targeting real hydrothermal alterations' zones within the study area. A great harmony between the detected alterations and our field observations was noticed (Figs. 15 and 16). Our field examination revealed different types of hydrothermal alteration minerals or already mined zones, manifesting the role of current research in deciphering the hydrothermal alteration pattern of the study area.
According to the geological setting (Fig. 1), El Rukham-Mueilha district is occupied by various lithologies, including suprastructural Neoproterozoic Pan-African nappe assemblage (dismembered masses of ophiolitic mélange and island arc association) and syn to post-orogenic intrusive rocks which represents the main geological units within the Barramiya-Mubarak sector 88 . These lithologies are heavily altered by the magmatic fluids in the name of hematization (Fig. 16a), serpentinization (Fig. 16b,c), argillic (Fig. 16d), and propylitic alterations (Fig. 16e). These various types of alterations are linked with magmatic activity accompanied by the formation of ophiolites and island arc association, which were synchronously emplaced by multi-phase magmatism, which allowed the transfer of hydrothermal fluids and, consequently, the leaching of parent rocks.
The nappe assemblage is the most prevalent lithology and is made up of serpentinites (Fig. 16), a metagabbrodiorite complex, and metavolcanics embedded within a matrix of low-grade metasedimentary rocks. Serpentinite represents the low-temperature metamorphic/hydrous alteration product of igneous ultramafic rocks, while the metavolcanics are mainly andesitic in composition. This greatly matches our results that highlighted serpentinite rocks within the study area as hydroxyl-bearing alteration zones represented mainly by Gabal Um salim and Gabal Um salatit (Fig. 15). The island arc metavolcanics are mainly made of andesite and are affected by hydrothermal alterations, especially in the southwestern part of the study area around wadi Dungash (Fig. 15). Petrographically phenocrysts of heavily altered plagioclase and hornblende that have completely altered into epidote, iron oxides, or chlorite and are embedded in a fine-grained groundmass of plagioclase, hornblende, quartz, and mica, giving it a porphyritic texture (Fig. 17a-e).
This sequence is later intruded by multi-stage magmatism, including syn-orogenic and post-orogenic granites. These syn-orogenic granitoids are mostly exposed in Gabal El-Rukham, as opposed to the post-orogenic granite that is intermittently encountered in the research region (i.e., Mueilha granite). The syn-orogenic granites are dioritic in composition, including plagioclase, hornblende, quartz, and epidote, which represent with iron oxides the alteration product of hornblende (Fig. 17d,f).

Discussion
Based on the previous hydrothermal alteration inspection using multi-sensor remote sensing data and airborne geophysical data, the main structural context of the study area, and previous studies 10,24,49,50,89 , the whole alteration pattern resulting from the four adopted datasets highlights their eligibility for detecting lithological units enriched with different types of hydrothermal alterations and their associated mineral deposits. After checking the eligibility of our results and coinciding with the main target of the current research, picking out the best sensor/s (based on their ability to delineate real alteration zones) to save time and effort, narrow the zones to be further exploited, and recommend them for further geological explorations is a critical task.  (Table 3) of the sensor results to select the best results. The study strongly recommends blending at least two optical sensors for confirming the altered zones and narrowing the area to be further investigated. For example, hydroxyl-bearing minerals are disclosed in the range of 2.10-2.28 μm. If this range is covered by different wavelengths and bandwidths (e.g., Landsat 8 b7 187 nm, Sentinel 2 b12 180 nm, Sentinel 2 b12 180 nm), it is far better to confirm the presence of these minerals. This is clearly shown by the results   Fig. 12, where the overestimated area is larger than the confirmed pixels (resulting from the intersection of the four resultant patterns of the utilized sensors). Areal differences in the estimated zones directly result from changes in band designations and specifications for each sensor (e.g., the covered wavelength and bandwidth). Even for the same sensor, multi-temporal data could result in various anomalies 7 , thus a combination of two optical sensors is strongly recommended for further geological investigations to obtain focused real anomalies. It should be emphasized that the overestimated areas may have real altered pixels. However, the presence of false alteration zones has a higher probability compared to zones confirmed by the four datasets (altered pixels are often over or among the mineralized spots compared to overestimated areas, which have many distal pixels as shown in Figs. 12, 13, and 14). Thus, starting field investigations with these narrow zones is recommended to save the cost of further geophysical and geochemical analysis. Moreover, these zones could guide us to some parts of the overestimated areas through structural elements or any other evidence during field investigations.
In the current research, ASTER data potentiality in alteration mapping is unparalleled with the other three sensors. Without a doubt, ASTER decommissioning from 2009 has a negative effect on the geological remote sensing community. As a solution for confirming ASTER data anomalies and selecting another appropriate sensor to be combined with ASTER in further investigations, anomaly overlay analysis of the other 3 datasets was executed with ASTER data. It should be emphasized that the similarity between ASTER and Sentinel 2 results is eye-striking (Figs. 3 and 5) and strongly supports the geologic model of the hydrothermal mineral deposits distributed within the study area, meaning that Sentinel 2 data could be applied in performing a vigorous geological mapping. Additional spatial investigations of the results with the approximate locations of the mineral deposits (Fig. 15) corroborated the efficacy of fusing Sentinel 2 data (that has shown its efficacy in highlighting comparable zones affected by hydrothermal alteration) with ASTER data using an anomaly-overlaying selection technique. In this way, Sentinel 2 data could revive ASTER data and help solve the time gap issue that often arises with ASTER's current investigations and reopen the door for applying the powerful SWIR spectral ranges in geological explorations. Moreover, magnetic and radiometric data implications clearly demonstrated how well the detected alteration anomalies of ASTER and Sentinel 2 are structurally controlled and mostly coincided with higher magnetic and K/eTh ratio anomalies.
It should be emphasized that the four datasets are eligible for comprehensive lithological and hydrothermal alteration mapping. However, In the current research, the favorability of ASTER and Sentinel 2 over L8 and ALI was assigned based mainly on the reality of the detected hydrothermal alteration zones (after comparison with airborne geophysical data and field investigations) which were attributed to the preferences of spatial and spectral resolutions as highlighted in Table 1. The latter highlighted that sentinel 2 data has 10 spectral ranges with a fine

Conclusions
The current research integrated multi-sensor remote sensing data, airbornemagnetic and radiometric data, field observations, and petrographical examinations to decipher the hydrothermal alteration pattern of the study area, highlighting the efficiency of these multispectral datasets, and recommending the best for the geological community. Cross-linking results of ASTER, Sentinel 2, Landsat 8, Advanced Land Imager, magnetic, and radiometric data in geological and hydrothermal alteration mapping revealed the following.   www.nature.com/scientificreports/ 1. The four datasets are eligible for comprehensive lithological and hydrothermal alteration mapping with acceptable areal differences attributed mainly to the powerfulness of the utilized wavelength and the desired target. 2. By comparing the results of at least two sensors, real alteration anomalies could be confirmed. This could narrow the area that may need to be investigated further with expensive geophysical and geochemical tests. Fieldwork should be commenced in these highlighted areas, which could lead to other anomalies that have not been found yet through other means, such as structural elements. 3. For a more detailed spectral analysis, ASTER data seems to be the best choice among the utilized sensors; however, a spatial analysis of the current study strongly recommends blending Sentinel 2 data results with ASTER findings to, confirm the altered zones, and narrow the area to be surveyed. 4. Band ratios, DPCA, and CEM techniques are extremely useful for deciphering the alteration patterns from the implemented datasets. Their results were pretty close to matching the structure of the area they looked at and how the mineralized spots were spread out. The current research strongly recommends integrating Sentinel 2 and ASTER data for deciphering further hydrothermal alteration patterns. Additionally, the placer deposits within the study area are of great interest due to the dominance of several alteration minerals in them. 5. Using the recommended integrated Sentinel 2 and ASTER data, the alteration zones were found to be in line with magnetic and radiometrically changed areas (with high magnetic and high K/eTh ratio anomalies) that are related to metalliferous minerals (like magnetite) and are located along the defined shear zone in the study area. Depth estimation analysis revealed that most of the anomalies have shallow sources (0-600 m); however, deep (up to 1500 m) structurally related anomalies also exist. 6. The present investigation unravels the hydrothermal alteration patterns (dominated mainly by OH-bearing minerals and iron oxides) and structural framework of the surveyed region by utilizing diverse data sources and scales, including satellite remote sensing images, airborne geophysical data, fieldwork, and petrographic examinations. Our research could enhance the likelihood of identifying genuine alteration zones for subsequent feasibility studies and exploitation. Huge argillic-ferric alteration zone, several tens of meters across, in serpentinites leading to brick red deposits that are mined by the random (artesian) miners with rather low 1-2 ppm gold according to their estimation. (e) Intensive random mining in the alluvial and highly fractured loose (greenish propylitic altered) serpentinites was done by the random (artesian) miners using the loaders and guided by the metal detectors. (All the participants identified in these field photographs are authors of the current research. These photos are our own and we agreed to publish them.)

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.